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Abstract 

Two-scalar theories at high temperature exhibit a rich spectrum of possible 
critical behaviour, with a second or first order phase transition. In the vicinity of 
the critical temperature one can observe critical exponents, tricritical points and 
crossover behaviour. None of these phenomena are visible to high temperature 
perturbation theory. 
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1. Introduction 



Scalar field theories have been the prototype for investigations concerning the question of 
symmetry restoration at high temperature. Following the original argument of Kirzhnits 
and Linde 0], the 0(iV)-symmetric scalar theory was considered in subsequent studies 
of the problem || [| |4|]. The framework in which these studies were carried out is the 
perturbative evaluation of the effective potential [||] and its generalization for non-zero 
temperature. Even though the restoration of the spontaneously broken symmetry was 
qualitatively demonstrated, the investigation of the details of the phase transition was 
not possible, due to infrared divergences rendering the perturbative approach unreliable 
near the critical temperature 0, f|]. These divergences originate in the absence of an 
infrared cutoff in higher loop contributions when the temperature dependent mass of the 
scalar fluctuations approaches zero near the critical temperature. An amelioration of 
the situation was achieved through the summation of an infinite subclass of perturbative 
contributions (the "daisy" graphs) 0. Indeed, these contributions become dominant 
for large iV and a quantitative description of the phase transition can be obtained in 
this limit. However, the physical picture remained unclear for small, physically relevant 
values of N, for which even the order of the transition was not established. The question 



was resolved || through the method of the effective average action @ - Unj? which relies 
on the renormalization group approach. The phase transition was shown to be second 
order for all values of N. The quantitative behaviour near the critical temperature was 
studied in detail and the critical system was found to have an effectively three-dimensional 
character. Its behaviour can be characterized by critical exponents || [TOf, in agreement 
with known results from three-dimensional field theory. The picture was verified through 
an independent analysis in the large N limit [|ll|], with use of other non-perturbative 
methods, such as the solution of the Schwinger- Dyson equations. A summary of the 
results can be found in ref. [fl2 1 . 

In a cosmological context first order phase transitions are more spectacular than sec- 
ond order transitions due to the departure from thermal equilibrium. One would like 
to have a prototype model for a first order transition, for which the methods of high 
temperature field theory can be tested, similarly to the O(N) scalar model for second 
order transitions. In statistical physics it is well known that scalar models with more 
than one field and discrete symmetries instead of maximal O(N) symmetry exhibit a rich 
spectrum of critical behaviour, including first and second order transitions and tricritical 
behaviour in between. Since high temperature field theories are in close correspondence 
to (three-dimensional) statistical models, it seems natural to investigate such models also 
as prototypes for first order transitions in high temperature field theory. In this paper 
we apply the method of the effective average action to the study of the high temperature 
phase transitions in theories with two real scalar fields. The symmetry is not 0(2), but 
rather a discrete symmetry. This model can serve as a prototype for a first order phase 
transition in field theories. It can be easily generalized to models in which each scalar 
field is an iV-component vector. 

We are interested in the phenomenon of spontaneous symmetry breaking and sym- 
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metry restoration at high temperature. For sufficiently low temperature our two scalar 
theory models the Higgs mechanism in gauge theories, through which the expectation 
value of a scalar field results in a mass term for gauge fields. f] Perturbative arguments 
predict a first order phase transition for this case However, the reliability of such 
predictions is questionable when the transition becomes weakly first order, due to infrared 
divergences similar to the ones plaguing the study of the 0(iV)-symmetric scalar theory 



14| . The approximate vanishing of some mass near the critical temperature results in 
the absence of an infrared cutoff in higher loop contributions of perturbation theory. In 
two-scalar theories this is connected with the fact that one of the fields gets its mass (or 
part of it) through the expectation value of the other. 

The present work obtains control over these infrared problems. Depending on the 
couplings of the model we find that the phase transition is either first or second order. 
For a sufficiently strongly first order transition high temperature perturbation theory 
may give realistic results for the two-scalar model. We concentrate here on the more 
problematic regions of a second order transition, a weakly first order transition and the 
tricritical behaviour at the separation of the two regimes. For the corresponding values 
of the couplings high temperature perturbation theory fails near the critical temperature. 
As a byproduct, our results can be used in order to establish in which region of parameter 
space perturbation theory gives a reasonable approximation for the description of the 
phase transition. 

Our results are directly relevant for two specific classes of scenaria in the cosmological 
context: The first class concerns multi-Higgs-scalar extensions of the standard model at 
non-zero temperature. The prediction of perturbation theory for a first order electroweak 
phase transition ||13|| , combined with the existence of baryon number violating processes 
at non-zero temperature within the standard model |18|], has generated much interest in 
the possibility of creating the baryon asymmetry of the universe during the electroweak 
phase transition ||19|| . It is not clear, however, if the phase transition in the pure standard 
model is sufficiently strongly first order and if there is sufficient CP violation in order to 
create an asymmetry of reasonable size. This has led to the study of multi-Higgs-scalar 
extensions of the standard model, in which the additional scalar fields can be used to 
make the phase transition more strongly first order or to enhance the sources of CP 
violation in the model. Also supersymmetric extensions of the standard model contain 
two scalar doublets. It is not clear whether the perturbative methods used in ref. |^D| for 
the calculation of the scalar field contributions to the effective potential are reliable for 
such two-scalar models, if we take into account the "warning" from the study of the O(N)- 

2 For sufficiently small gauge coupling the present investigation and its generalization to the case where 
each scalar field has N components gives a reasonable approximation to the gauged models even in the 
vicinity of the critical temperature. However, the determination of the meaning of "sufficiently small" 
needs a detailed investigation of high temperature gauge theories. 

3 We should point out that, for a non-abelian Higgs model, the situation is much more involved than 
the perturbative results indicate, due to the presence of a confining regime in the symmetric phase of the 
model. As this work deals only with scalar fields, such a complication does not arise. For a discussion of 



gauge theories in the context of the effective average action approach see refs. |15, 16l uj\ 
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symmetric theory. Our work gives a reliable estimate of the effect of these contributions 
on the nature of the transition. 

The second class of scenaria concerns multi-scalar models of inflation ||21J| . I n most 



such studies some classical potential is employed, which may bear no resemblance to the 
effective potential. Thermal effects are often ignored except for the temperature depen- 
dence of the mass term. If inflation is initiated by a high temperature phase transition 
our formalism sets the framework for the proper study of the problem. 

We consider a theory of two real scalar fields Xa{ a — 1 ; 2), invariant under the discrete 
symmetries (xi <-» — Xi> X2 ~X2, Xi X2), which we denote by (1 «-> —1, 2 <-> —2, 1 «-> 
2) for brevity. The symmetry group is Z4 x Z 2 , consisting of 90° rotations in the (xi, X2) 
plane and a reflection on one of the axes. The classical potential can be written as 

v( X i, X2) =\m\xl + xl) + \kx\ + xi) + \9xlxl 

= \m 2 {xl + Xl) + \Mxi + Xl? + \x\xlxl (1.1) 

with 

x=j-l- (1.2) 

For V to be bounded from below we require A > 0, x > —2. 

For m 2 > the classical theory is in the symmetric regime (which we denote by S) 
with the minimum of the classical potential at the origin. For m 2 < the theory is in 
the spontaneously broken regime and we distinguish two possibilities consistent with the 
symmetry: 

I) For x < four degenerate minima of the potential are located between the two axes at 



/ 2m 2 

Xio = ±X20 = ±\- T —- (1.3) 
V X + 9 

We denote this regime by M. 

II) For x > the four minima of the potential are located on the axes at 



/ 2m 2 

Xio = ±y j-, X2o = 0, (1.4) 

or similarly with Xw an d X20 interchanged. We denote this regime by AX. 

The regimes M and AX are closely related. A redefinition of the fields according to 

Xi = -j={Xi + X2) 

X2 = -^[Xi-X2) (1.5) 

results in a rotation of the axes by 45°, thus transforming the AX into the M regime. The 
couplings of the redefined theory are related to the old ones according to 

A = Xfl + f 
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* = -TT^ (1 - 6) 

There are three characteristic values of x: 

a) For x = the symmetry of the theory is increased to 0(2). 

b) For x = — 1 the theory decomposes into two disconnected Z 2 -symmetric models for \i 
and X2 separately. 

c) Similarly, for x = 2 (and therefore x — — 1, according to eqs. (|1.6|) ) the theory decom- 
poses into two disconnected Z 2 -symmetric models for xi an d Xi- 

The above symmetries are expected to be preserved after the quantum or thermal cor- 
rections have been taken into account. This means that any renormalization group flow 
of the couplings that starts on the surfaces x — 0, —1,2 in parameter space cannot take 
the system out of them. As a result the parameter space of the theory is divided into the 
four regions: x>2, 2>x>0,0>x> — 1, x < — 1, which are not connected by the 
renormalization group flow of the couplings. The phase transitions for the theories which 



correspond to x = 0, —1, 2 have been discussed in detail in refs. || |10| . They are second 
order transitions governed by effectively three-dimensional fixed points. In our model 
these fixed points exist on surfaces separating the parameter space into disconnected 
regions. We shall demonstrate all the above points in the following sections. 

We should point out that this model was discussed in ref. [E^] through use of finite 
temperature perturbation theory. No part of the rich structure of critical behaviour that 
we shall describe in the following sections was observed. The universal, effectively three- 
dimensional behaviour of the system near the critical temperature is common for statistical 
systems and three-dimensional field theories which belong to the same universality class. 
(The statistical systems are characterized as two-component spin systems with cubic 
anisotropy.) As a consequence, various aspects of the problem have been investigated in 
refs. |23f - [^(| (and references therein) through other methods. Our results are in very 



good agreement with all these studies. Similar models have been considered in ref. [27 



The outline of our procedure is the following: We make use of the effective average 
action Tk, which results from the effective integration of quantum and thermal fluctuations 
with characteristic momenta q 2 > k 2 . It contains all the information on the generalized 
couplings of the theory and their dependence on the scale k. For k of the order of some 
ultraviolet cutoff A the effective average action is equal to the classical (bare) action 
(no integration of fluctuations takes place). For k = 0, is equal to the effective 
action (all fluctuations are integrated). The dependence of T k on the scale k is given 
by an exact non-perturbative renormalization group equation, which can be expressed 
as evolution equations for the running couplings of the theory. These equations can be 
solved within some appropriate approximation scheme, with the classical couplings as 
initial conditions for k = A. The renormalized couplings of the theory are obtained 
for k = 0. The calculation can be performed for zero and non-zero temperature. The 
gradual incorporation of the effects of quantum and thermal fluctuations into the running 
couplings is the essential element which resolves the problem of infrared divergences that 
invalidates perturbative schemes. The basic formalism of the effective average action is 
summarized in section 2. We expect that the running of the couplings for k much smaller 
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than the temperature has an effectively three-dimensional character. The reason is that 
the effective dimensionality is reduced, when the characteristic length scale 1/k of the 
"coarse grained" system is much larger than the periodicity 1/T in the imaginary time 
direction set by the temperature. This is expected to be important near the critical 
temperature, where no infrared cutoffs (such as masses) other than k exist. As a result, 
the fixed point structure of the three-dimensional theory determines the behaviour of 
the critical system. For this reason we present in section 3 a qualitative study of the 
three-dimensional theory and its fixed points, based on a crude approximation scheme for 
the solution of the exact renormalization group equation. In section 4, we develop more 
elaborate (and, therefore, more accurate) approximation schemes. They are generalized 
for non-zero temperature in section 5. These tools are put into work in sections 6 and 7: 
The evolution of the running couplings is calculated, starting with the classical theory at 
scales k = A ^> T and finishing at k — 0, where the renormalized theory is obtained. In 
section 7 we explicitly demonstrate how the evolution of the running couplings becomes 
effectively three-dimensional for k < T. In section 8 we calculate the critical temperature 
for the phase transition. In sections 9-11 we discuss the details of this transition. We 
observe a rich spectrum of critical behaviour with critical exponents, crossover phenomena, 
tricritical points etc. None of these are visible within perturbation theory. Our conclusions 
are given in section 12. 



2. The evolution equation for the effective average 
potential 

We consider a theory of two real scalar fields Xa{ a — 1,2), in d- dimensional Euclidean 
space, with an action S[x] invariant under the (1 <-> —1,2 «-> —2, 1^2) symmetry. We 
specify the action together with some ultraviolet cutoff A, so that the theory is properly 
regulated. We add to the kinetic term an infrared regulating piece || 

AS = \ I -^-Ml)x:(q)x a (q), (2-1) 

where x a (<?) are the Fourier modes of the scalar fields. The function is employed in 
order to prevent the propagation of modes with characteristic momenta q 2 < k 2 . This 
can be achieved, for example, by the choice 

with 

flid) = exp • (2.3) 

We point out that there are many alternative choices for Rk{q), some of which were used 



in refs. @ - []12| . The physical results which are obtained when the cutoff is removed 
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are scheme independent. The choice of eqs. ( |2.2| ),( PTB| ) is the most natural one || and 
convenient for numerical calculations. For a massless field the inverse propagator derived 
from the action S + AS has a minimum ~ k 2 . The modes with q 2 ^> k 2 are unaffected 
by the infrared cutoff, while the low frequency modes with q 2 <C k 2 are cut off, as R k acts 
like a mass term 

lim R k (q) = k 2 . (2.4) 

We subsequently introduce sources and define the generating functional for the connected 
Green functions for the action S + AS. Through a Legendre transformation we obtain 
the generating functional for the 1PI Green functions Ffc[0 a ], where <p a is the expectation 
value of the field \ a m the presence of sources. The use of the modified propagator for the 
calculation of r\ results in the effective integration of only the fluctuations with q 2 > k 2 . 
Finally, the effective average action is obtained by removing the infrared cutoff 

t~i tJM S r, al 1 f d d q 



I ^M^WiQ)- (2-5) 



J 2 J (2tt) 

For k equal to the ultraviolet cutoff A, becomes equal to the classical action S (no 
effective integration of modes takes place), while for k — > it tends towards the effective 
action T (all the modes are included) which is the generating functional of the 1PI Green 
functions computed from S (without infrared cutoff). For intermediate values of k the 
effective average action realizes the concept of a coarse grained effective action in the 
sense of ref . |28[ . 



The interpolation of between the classical and the effective action makes it a very 
useful field theoretical tool. The means for practical calculations is provided by an exact 
flow equation Q which describes the response of the effective average action to variations 
of the infrared cutoff (t = ln(fc/A)) @ 

|r fc M = ^Tr{(rl 2) M + R k y l jR k }. (2.6) 

functional derivative of tl 
to cf) a . For real fields it reads in momentum space 



(2) 

Here T k is the second functional derivative of the effective average action with respect 



with 

r(-q) = €(q)- (2-8) 

The non-perturbative flow equation has the form of an one-loop expression involving 
the exact inverse propagator T k together with an infrared cutoff provided by R k . No 
contributions from higher loops appear in this exact equation. 



See ref. p9] for other versions of exact renormalization group equations. 
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For the solution of eq. ( |2.6p one has to develop an efficient truncation scheme. The 
form of the effective average action is constrained by the (1 «-> —1,2 <-> — 2,1 <-> 2) symme- 
try. However, there is still an infinite number of invariants to be considered. Throughout 
this paper we shall work with an approximation which neglects the effects of wave function 
renormalization. We shall, therefore, keep only a classical kinetic term in the effective 
average action 

T k = J d d x{u k ( Pl ,p 2 ) + ^dy a d^ a }, (2-9) 

and neglect all invariants which involve more derivatives of the fields. We have used the 
definition p\ = ^<p\ and similarly for p 2 . The justification for our approximation lies in 
the smallness of the anomalous dimension, which is expected to be r\ ~ 0.03 — 0.04 for the 
three-dimensional theory. We estimate the corrections arising from the proper inclusion 
of wave function renormalization effects to be of the same order as rj (a few %). An 
improved treatment will be given elsewhere In order to obtain an evolution equation 



for Uk from eq. ( |2.(j| ), we have to expand around a constant field configuration (so that 
the derivative terms in the parametrization ( |2.9| ) do not contribute to the l.h.s. of eq. 
Eq. O then gives § - []TD| 



lu k ( Pl , P2 ) = I J ^- ( \ + p( 2 * ) (2-10) 



(2n) d \P(q 2 ) + M? P{q 

P{q 2 ) results from the combination of the classical kinetic contribution q 2 and the regu- 
lating term R k into an effective inverse propagator (for massless fields) 

p(q2) = q2 + Rk = T^m Y (2 - n > 

with fl{q) given by eq. ( |2.3| ). For q 2 ^> k 2 the inverse "average" propagator P(q) 
approaches the standard inverse propagator q 2 exponentially fast, whereas for q 2 <C k 2 
the infrared cutoff prevents the propagation. Mf 2 are the eigenvalues of the mass matrix 
at the point (pi, p 2 ) 

Mh(pup 2 ) = l\Ux + U 2 + 2U llPl + 2U 22 p 2 



2 



± 



(Ut -U 2 + 2U nPl - 2U 22 p 2 ) 2 + lW 2 12 p lP2 1 , (2.12) 



and we have introduced the notation Ui = U\ 2 = ^ 2j i k etc. 

1 api ' ±z op\op2 



Eq. ( |2.10| ) is the master equation for our investigation. It is a non-linear partial 



differential equation for three independent variables {t,p\,p 2 ). Since it is difficult to solve 
it exactly we again resort to some approximation scheme. We first introduce a Taylor 
expansion of U k (pi, p 2 ) around its minimum. This turns eq. ( |2.10|) into an infinite system 



of ordinary (coupled) differential equations for the fc-dependence of the minimum and the 
derivatives of the effective average potential, with independent variable t = ln(fc/A). We 



7 



solve this system approximately by truncating at a finite number of derivatives. This 
approach has been used in the past for the study of the 0(iV)-symmetric scalar theory. 
It has provided a full, detailed picture of the high temperature phase transition for this 
theory || , [|lOj - Jl2| , with accurate determination (at the few % level) of such non-trivial 
quantities as the critical exponents [fLOfl . An estimate |30| of the residual errors for high 
level truncations indicates that they are smaller than the uncertainties introduced by the 
imprecise treatment of the wave function renormalization effects. For this work we shall 
use the lowest level truncation which keeps only the second derivatives of the potential 
Un, U 22 , U\2- This will be sufficient for a reliable determination of the phase diagram and 
a crude estimate of universal quantities such as critical exponents and crossover curves. 
For an improved treatment see ref. p3 |, and for a discussion which takes into account 
the next level in the truncation for Uk and the first corrections arising from wave function 
renormalization see ref. p3 1 . 



3. The phase structure of the three-dimensional the- 
ory 

Before performing a more detailed analysis we would like to gain some understanding of 
the phase structure of the theory. As we have already mentioned in the introduction, 
the behaviour of the four- dimensional theory near a high temperature second order phase 
transition is expected to have a three-dimensional character. The reason for this is the 
divergence of the correlation length for the fluctuations of the system (the mass of some 
field goes to zero). As a result, the characteristic length scale for the critical system is 
much larger than the periodicity in the imaginary time direction due to temperature (for 
details see the following sections). | This leads to dimensional reduction and the critical 
system has effectively three-dimensional behaviour. For this reason we are interested 
in the phase structure of the three-dimensional theory. More specifically we want to 
investigate the possible existence of fixed points which govern the dynamics of second 
order phase transitions. 

For the purpose of this section we parametrize the potential by its derivatives at the 
origin (S regime) 

m\k) = I7i(0) = U 2 (0) 
Mk) = U n (0) = U 22 (0) 

g(k) = u 12 (o) 

m = f| - 1. (3.1) 

The equality of Ui, U 2 and Un, U 22 is imposed by the (1 <-> —1,2 <-> —2,1 «-> 2) 
symmetry of the theory. For the potential to be bounded we also require x > —2. For 

5 If the phase transition is strongly first order this need not be true, because the mass of the fluctuations 
does not go to zero and the correlation length does not diverge near the transition. 
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a rough estimate the three-dimensional couplings are related to the effective couplings 
of the four- dimensional theory at high temperature by A(27rT) = A4T, g{2irT) = g 4 T ', 
m 2 (2irT) = m\ + cT 2 , with appropriate c (for details see sections 6 and 7). Evolution 
equations for the above parameters are obtained by taking derivatives of eq. (|2.10| ) with 
respect to pi j2 . The truncated flow equations read for d = 3 

— = v 3 k(4 + x)XL 3 1 (m 2 ) (3.2) 

^ = -^AT^IO + 2x + x 2 )X 2 L 3 2 (m 2 ) (3.3) 

^ = v 3 k- 1 (x + l)x(x-2)XL 3 2 (m 2 ), (3.4) 

with v 3 = 1/87T 2 . The threshold functions L^(w) suppress the contributions of massive 
modes to the evolution equations. They are studied in detail in the following sections and 
in appendix A. It is convenient to define the dimensionless couplings 



m 2 {k) 
X(k) 



m 2 {k) 
k 2 

A(fc) 
k 



g{k) = (3.5) 

In terms of these quantities the evolution equations have a scale invariant form, in the 
sense that the r.h.s. does not explicitly involve a dependence on k 

— = -2m 2 + v 3 {4 + x)\Ll(m 2 ) (3.6) 

^ = -X-v 3 (10 + 2x + x 2 )X 2 L 3 2 (m 2 ) (3.7) 

^ = v 3 (x + l)x(x-2)XLl(m 2 ). (3.8) 

We are interested in the fixed points of the last set of equations. For any x, there 
is an ultraviolet attractive Gaussian fixed point with m 2 = X = 0. There are also three 



fixed points with at least one infrared attractive direction [p3] - 0]. They all appear 
for m 2 < 0, A > 0. (The exact values are not important since the discussion in this 
section is only qualitatively correct.) For their identification we use their standard names 
in statistical physics [^3[ 

a) The Heisenberg fixed point has x = and corresponds to a theory with symmetry 
increased to 0(2), as we have discussed in the introduction. 

b) The Ising fixed point has x = — 1 and corresponds to two disconnected ^-symmetric 
theories. 
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c) The Cubic fixed point has x = 2 and corresponds to two disconnected theories, if the 
fields are redefined according to eqs. fll.5|) . 

All these points are infrared unstable in the m 2 direction and are located on a critical 
surface m 2 cr = ml r (\,x) < 0. Solutions of the evolution equations which start above the 
critical surface, with m 2 > m 2 r1 flow towards the region of positive m 2 for t — > — oo, 
and correspond to theories in the symmetric phase. Solutions with m 2 < m 2 r flow deep 
into the region of negative m 2 and correspond to theories in the phase with spontaneous 
symmetry breaking. 

The relative stability of the fixed points on the critical surface determines which one 
governs the dynamics of the phase transition very close to the critical temperature. For 
a first simple investigation of the relative stability in the (A, x) directions we fix m 2 
to an arbitrary value (we choose m 2 = for convenience) and solve eqs. ( |3.7| ), ( |3.8|) 
numerically. The results are presented in fig. 1. We observe that all three fixed points 
are attractive in the A direction. However, the Ising and Cubic fixed points are repulsive 
in the x direction, while the Heisenberg fixed point is totally attractive. We observe four 
disconnected regions: 

a) 2 > x > : The trajectories flow away from the Cubic and towards the Heisenberg 
fixed point. 

b) > x > — 1 : The trajectories flow away from the Ising and towards the Heisenberg 
fixed point. 

c) x > 2 : The trajectories flow away from the Cubic fixed point and into a region 
of large x and small A. Eventually A turns negative at a finite value of k. (This can 
be verified through the explicit solution of eqs. (|3.7| ), ( |3.8|) in this region.) At this 
point an instability arises, as the potential seems not to be bounded from below. Our 
treatment is not sufficient for a detailed investigation of the nature of this instability, since 
our truncation scheme is very crude. A detailed discussion is given in ref. ||33|| , where 
improved truncations are employed. It is shown that the instability is not real since the 
higher derivatives of the potential remain positive. The change of sign for A corresponds 
to the disappearance of a false vacuum of the theory and results in a first order phase 
transition. We shall return to this point in the following sections. 

d) x < — 1 : The trajectories flow away from the Ising fixed point and cross the line x = —2 
at a finite k. This again implies the presence of an instability whose true nature is related 
to the disappearance of a false vacuum. The model exhibits a first order transition also 
for x < — 1. 

Flows that start on the lines x = 0,-1,2 remain on these lines. No trajectories exist 
which connect the four regions x>2, 2>a;>0, 0>x> —1, x < —1. All this is in 
agreement with the discussion at the end of the introduction. 

The diagram of fig. 1 determines the phase structure of the theory when the behaviour 
of the system becomes effectively three-dimensional (i.e. close to the critical temperature). 
For parameters in the regions 2 > x > 0, > x > -1 we expect second order phase 
transitions, with critical dynamics governed by the three fixed points. These two regions 
can be mapped onto each other through a redefinition of the fields according to eqs. (|1.5|) , 
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Ql.6| ). This indicates that the Ising and Cubic fixed point should lead to identical universal 
behaviour (and therefore to identical universal quantities, such as critical exponents). Very 
close to the critical temperature we expect the Heisenberg fixed point to dominate the 
dynamics, but the other two can be relevant if the initial values of the running parameters 
are sufficiently close to them. In the parameter regions x > 2, x < — 1 we expect first 
order phase transitions. In the following sections we shall verify the above conclusions 
with improved quantitative accuracy. 



4. Truncations of the evolution equation 

We proceed now to a more detailed study of the evolution equation and its truncations. As 
we have discussed at the end of section 2, we parametrize the effective average potential by 
its minimum and its derivatives at the minimum. For this work we shall use a truncation 
that preserves up to second derivatives of the potential. 

In the symmetric regime (which we denote by S) the minimum of the potential is at 
Pio(^) = P2o(k) = and we use the definitions of eqs. ( |3.1| ). The evolution equations in 
arbitrary dimension d preserve automatically the symmetry. They read 

' //7r> v d k d ~ 2 (A + x)XL d (m 2 ) (4.1) 



dt 



^ = - Vd k d -\10 + 2x + x 2 )\ 2 L d (rh 2 ) (4.2) 
^ = v d k d -\x + l)x(x - 2)\L d (m 2 ), (4.3) 



with the dimensionless integrals L d (w) given by 

'd\ r dP 



L d n (w) = - nk 2n - d n-i T 1 1 \J d d q^-(P + w)~^ 



- nk 2n ~ d I" dxxi-^{P + wY {n+l) . (4.4) 



Here P is given by eq. (|2.11 ), and 



vj 1 = 2^vrir (jj . (4.5) 

In the spontaneously broken regime there are two possibilities consistent with the sym- 
metry: 

I) In the M regime the minimum of the potential is located symmetrically between 
the p-axes at pw(k) = p2o(k) = ~p (k). We define the couplings 

X(k) = U n ( Po ) = U 22 (p ) > 

g( k ) = U 12 { P0 ) 

xik) = f| - 1. (4.6) 
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The requirement that the point (|po, \po) is the minimum of the potential imposes x < 0, 
while the potential is bounded at infinity for x > —2. For x = the symmetry of the 
theory is increased to 0(2) and the potential develops a series of degenerate minima 
along the circle pio + P20 = Po- For x — — 1 the theory decomposes into two disconnected 
^2(01,2 ^ — 0i i 2)-symmetric models. The mass eigenvalues are given by Ml = (2 + x)Xpo, 
M| = — x\pQ. The evolution equation for the minimum po(k) is obtained by considering 
the total t-derivative of the conditions fnflpo = f^flpo = |, |- Again, the truncated 
evolution equations automatically preserve the symmetry and read 

^ 2 {3Lf((2 + x)Xp ) + |^Lf(-xApo)} (4.7) 



^= -^- 2 {3L?((2 + x)Ap ) + | T |- 

§ = - ^- 4 ^i±l{L?((2 + x)Ap ) - Li(-xXp )} 
at po 1 + x L J 

2 / ™\ 2 



fix 




1 + ^ L 2 d ((2 + x)Ap ) + (l - I) L d 2 (-x\p )j (4.8) 
+ {Lf((2 + x)Ap ) - Lf(-xApo)} 



1 + ^V^((2 + x)Ap ) + (l - ^ 2 ^(-xApo)l (4.9) 



2) zvv ' ' ' \ 2 



For x = the above evolution equations reproduce the equations of the 0(2)-symmetric 
theory, while for x = —1 those of the ^-symmetric one (compare with refs. || p~0|| ) . 



II) In the regime which we denote by AX, two degenerate minima of the potential exist 
on each one of the p-axes. Without loss of generality we concentrate on the minimum 
at pio(k) = po(k),p2o(k) = 0. At the level of truncations that we are considering, the 
remaining parameters of the theory are conveniently defined as 

X(k) — Uu(po) 
g{k) = U 12 (p ) 



x{k) 



X(k) 

m 2 2 (k) = U 2 (p ). (4.10) 

The symmetry demands that for the truncated potential 

m\(k) =x(k)X(k)p {k). (4.11) 

The requirement that the point (po, 0) is the minimum of the potential imposes x > 0. As 
before, for x = the symmetry of the theory is increased to 0(2). The mass eigenvalues 
are given by Mf = 2Xpo, Mf = a;Apo- At this point we encounter a difficulty. The deriva- 
tion of truncated evolution equations is hindered by the fact that the parametrization 
around a minimum located on one of the axes is asymmetric between the two fields. As 
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a result the symmetry (0i <-> ^2), is not maintained by the evolution equations at each 
level of the truncations. More specifically the flow equations for the couplings Un(p ), 
£^22 (po) are different. Also eq. (|4.11|) is not preserved by the evolution equation. This 
is not surprising, since these relations are not expected to hold for the exact potential 
without truncation. It is easy to see that they are altered as soon as third derivatives 
of the potential are included. This is in constrast with what happens in the M regime, 
where the formulation is symmetric between the two fields. For example, in the M regime 
the couplings Un, U22 are expected to remain equal at every level of truncations, and 
indeed this is guaranteed by the evolution equations. The above remarks indicate a nat- 
ural method of preserving the (1 — 1, 2 <-> —2, 1^2) symmetry in the AX regime. A 
redefinition of the fields according to 



)2 = T2 m 



(4.12) 



results in a rotation of the axes by 45°, thus transforming the AX into the M regime. The 
couplings of the redefined theory are related to the old ones according to 



Po 
X 

x 



Po 

111 



X 



1 + 2 
- 1 - ~ 2 



(4.13) 



The evolution equations for the new quantities have already been worked out and are 
given by eqs. 



( |4.9| ). By simply rewriting them in terms of the old quantities defined 
in eqs. ( [4.1 0|) we obtain 



df>o 
dt 
dX 



Vd 



-77 = - v d 



k d - 2 {3Li(2X Po ) + (1 + x)Li{x\p )} 
k d ^X 2 {9L d (2Xp ) + (1 + x) 2 L d 2 {xX Po )} 



(4.14) 
(4.15) 



dx 

H= Vdk 



d-2 



6 ' + i{Lt(2Xpo) 



Po 1 

+ v d k d - 4 xX{9L d 2 (2Xp 



-Lf(xApo)} 

[l + x) 2 L d 2 (xX Po )}. 



(4.16) 



For x = the above evolution equations reproduce the ones of the 0(2)-symmetric theory. 
Another special point is x = 2 for which the theory, when expressed in terms of the 
redefined fields 0i, <p 2 , decomposes into two disconnected ^2(^1,2 <-> — 0i j2 )-symmetric 
models. We should point out that eqs. (|4.14| )- (|4.16 ) could have been obtained by defining 
X = Un(po) and x = £/2(po)/^n(po)po (in agreement with eq. ( |4.11|) ) and inserting eqs. 
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( [4.1 Op on the r.h.s. of the flow equations. The advantage of the redefinition ( |4.12| ) is 



that it makes transparent how this apparently arbitrary choice of parameters preserves 
the original symmetry at this truncation level. 

In the following sections we shall use the evolution equations (|4J| )-(B, (Q-Q, 



( [4.14|) -( f4.16|) , for the S, M, AX regimes respectively, in order to obtain the renormalized 



theory in its various phases. 



5. The integrals L d n for zero and non-zero temperature 

The integrals L^(w), defined in eq. ( |4.4| ), have been discussed extensively in refs. [§, |H], 
pT| (for various shapes of the infrared regulating function Rk(q), for which eqs. (|2.2| ),( pT3|) 
are the most natural choice |J). For completeness, we give in appendix A a summary of 
the properties of L^{w). Also, as an example, we plot in fig. 2 the integrals £?(u>), L\{w). 
Their most interesting property, for our discussion, is that they fall off for large values 
of w/k 2 , following a power law. As a result they introduce threshold behaviour for the 
contributions of massive modes to the evolution equations. It is obvious from eqs. ( f4.1| )- 
(pTTD-(|4"Tg|), ( f4.14j )- Q4"TTS| ), for the S, M, AX regimes respectively, that the various 



contributions to the evolution equations involve L d n integrals with the mass eigenvalues 
as their arguments. When the running squared mass of a massive mode becomes much 
larger than the scale k 2 (at which the system is probed), these contributions vanish and the 
massive modes decouple. We evaluate the integrals L^w) numerically and use numerical 
fits for the solution of the evolution equations. 

In order to extend the formalism of the previous section to non-zero temperature 
we only need to recall that, in Euclidean formalism, non-zero temperature T results in 
periodic boundary conditions in the time direction (for bosonic fields), with periodicity 
1/T |32fl . This leads to a discrete spectrum for the zero component of the momentum go 



q ^2nmT, m = 0, ±1, ±2, ... (5.1) 

As a consequence the integration over go is replaced by a summation over the discrete 
spectrum 

^ -TZ[£&. (5.2) 



(21)" f J (27r)"- 1 ' 

With the above remarks in mind we can easily generalize our master equation fl2.10|) 
in order to take into account the temperature effects. For the temperature dependent 
effective average potential P2,T) we obtain 

l^.^^r^E/^-v ( ? J Sf+ ^)|« h (5.3) 

with the implicit replacement 
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in eqs. ( |2.2| ),( pTB| ) and ( J2.1 1|) for and P. Again, the usual temperature dependent 



effective potential f2|-f§] is obtained from £4(pi, p 2 , T) in the limit k — > 0. As before, we 
can parametrize Uk{pi,P2, T) in terms of its minimum and its derivatives at the minimum. 
The evolution equations are given by ( fOl )- (|4~3"D , (|4.7|) - (^9|) , fl4.14|) -( fI7L6l) , with the 
obvious generalizations 



Po(k) - 




X(k)- 


-A(fc,T) 


x(k) - 


-+x(k, T) 


m 2 (k) - 


^m 2 (fc,T) 



(5.5) 

The integrals for non-vanishing temperature read 

L d n (w,T) = -2nk 2n - d K-i +1 T (-) f d^q ^-{P + w)^ n+l \ (5.6) 

where the implicit replacement (|5 .4] ) is assumed in P. Their basic properties can be estab- 
lished analytically. For T -C k the summation over discrete values of m in the expression 
( |5.6| ) is equal to the integration over a continuous range of go up to exponentially small 
corrections. Therefore 

L d n {w,T) = L d n {w) for T < A;. (5.7) 

In the opposite limit T ^> k the summation over m is dominated by the m = contribu- 
tion. Terms with non-zero values of m are suppressed by ~ exp (—(mT/k) 2 ) . The leading 
contribution gives the the simple expression 

L d n {w,T) = V -^^-L d n -\w) forT»A;, (5.8) 

v d & 

with Vd defined in (|4.5[) . The two regions of T/k in which L d (w, T) is given by the equa- 
tions (|5.7|) , (|5.8| ) are connected by a small interval, in which the exponential corrections 
result in a more complicated dependence on w and T. 

The above conclusions are verified by a numerical calculation of Lf(w, T). In fig. 3 we 
plot L\{w, T)/L\{w) as a function of T/k, for various values of w/k 2 . We distinguish 
three regions: 

a) T/k < 6i\ This is the low temperature region where L\ 2 {w 1 T) are very well approx- 
imated by their zero temperature value. We take 6\ = 0.15 and use Lf 2 (w,0) in the 
evolution equations for k > T/Q\. 

b) Q\ < T/k < 02'. In the threshold region we perform a numerical fit of the curve corre- 
sponding to w = which we use for all values of w. This is a good approximation since 

6 Comparison with results presented in rcf. |6| shows that the form of these functions depends on 
the details of the infrared regulating function Rk(q)- However, the physical results, which are obtained 
when the cutoff is removed, are independent of the shape of the cutoff. This will be apparent in the next 
sections and is a verification of the scheme independence of our conclusions. 
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the relevant w/k 2 turns out to be small in this region (see next sections). 

c) T/k > 9 2 : We take 9 2 = 0.4. For the high temperature region we use for the numerical 

solution of the evolution equations 

Ll 2 (w,T)=4^Ll 2 (w). (5.9) 

The three dimensional character of the effective theory for modes with q 2 <C T 2 manifests 
itself in the appearance of the three dimensional momentum integrals. It acquires here a 
precise quantitative meaning. 

We have now developed the necessary formalism for the study of the four-dimensional 
zero and non-zero temperature theory. In the following two sections we study the evo- 
lution of the running parameters of the theory, which leads to the determination of the 
renormalized theory at zero and non-zero temperature. 

6. The running in the low temperature and threshold 
regions 

In section 4 we derived the zero temperature evolution equations for the parameters 
(masses, vacuum expectation values and couplings) of the truncated theory as a function 
of the scale k in the various regimes (S, M, AX). In section 5 we generalized the formalism 
in order to take into account non-zero temperature effects. The evolution equations can be 
solved for a given set of initial conditions, specified as the values of the running parameters 
at a scale equal to the ultraviolet cutoff of the theory (k = A). As we pointed out in section 
2, at this scale the effective average action is equal to the classical action. Therefore, the 
initial values for the parameters correspond to their classical (or bare) values. Also, 
the discussion in section 5 has shown that in the low temperature region {k > T/9i) 
there is no difference between the zero and non-zero temperature theory. As a result, we 
can define the theory in terms of the classical values of its parameters at k = A ^> T, 
independently of the temperature. The integration of the evolution equations gives the 
running couplings at lower scales. No temperature effects are observed in the evolution 
inside the low temperature region (k > T/9i). We can, therefore, use the values of the 
running couplings at k = T/9i for the definition of the theory, since they are in one-to- 
one correspondence with the classical couplings, independently of the temperature. This 
turns out to be the most convenient choice and we shall use it for the rest of the paper. 
The temperature starts to become important when the evolution enters the threshold 
region (T/9i > k > T/9 2 ). In the high temperature region {k < T/9 2 ) the evolution is 
effectively three-dimensional, as we discussed in section 5. Finally in the limit k — > the 
effective average action becomes the effective action, and the integration of the evolution 
equations gives the renormalized values for the couplings at various temperatures. All 
the information on the various phases of the theory is contained in these renormalized 
couplings and their temperature dependence. 
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We have seen in the introduction and in section 4 that the AX regime (x > 0) and 
the M regime (x < 0) can be mapped onto each other through a simple redefinition of 
the fields (see eqs. (|4.12|) , (|4.13| )). For this reason, the physical behaviour in the two 
regimes is the same. For example, the Cubic and Ising fixed points generate the same 
universal behaviour (characteristic of a Z 2 -symmetric scalar theory). For this reason, we 
shall concentrate on the region x > only. All the results can be easily extended to the 
region x < 0, through the transformations of eqs. ( [4.12|) , (|4.13|) . 

Since we are interested in symmetry restoration at non-zero temperature, we first 
consider the theory in the spontaneously broken regime. The evolution equations in 
the AX regime (which is the relevant one for x > 0) in four dimensions and non-zero 
temperature can be easily derived from eqs. ( |4.14|) - (|4.16| ) and read 



dpo 
dt 
dX 
~di 
dx 
~dt 



-^ 2 {3^(2Apo)ti(2Apo,T) + (l + x)^(xApo)ti(xAp ,T)} (6.1) 

- v 4 \ 2 {9L*(2\p )t 2 (2\p ,T) + (1 + x) 2 L%(x\p )t 2 (x\p , T)} (6.2) 

^4-^^{^(2Apo)ti(2Apo,T)-Lt(xApo)ti(xApo,T)} 

+ v 4 x\{9L 4 2 (2\p )t 2 (2\po, T) + (1 + x) 2 Lj{x\p )t 2 {x\p , T) j. (6.3) 



where = 1/327T 2 . We have not indicated explicitly the k and T dependence of the 
running parameters po(k,T), A(fc,T), x(k,T). They are defined at zero temperature 
according to eqs. (|4.10|) , and generalized for non-zero temperature according to eqs. 
750 • The functions ti 2 are defined as 



t n (w,T) 



L A n {w,T) 
Li(w) 



(6.4) 



with ti(w,T) plotted in fig. 3. 

At zero temperature one has t n (w,0) = 1 and the evolution equations have only one 
infrared attractive fixed point, the Gaussian one. In the limit of small A we shall neglect 
the slow logarithmic running of A, which is eventually stopped by the mass terms in 
the threshold functions L\ 2 . Similarly the running of x can also be neglected since it is 
suppressed by A/327T 2 . (For small A the difference of the two L\ functions in the first 
line of eq. ( |6.3| ) gives a contribution oc AL 2 (0).) For large A the evolution equations can 
be integrated numerically and the small resulting corrections can be reliably computed. 
This has been done in ref. |J for the 0(iV)-symmetric scalar theory. In this paper we 
concentrate on small couplings for which analytical expressions can be obtained. Eq. 
( |6.1|) can be integrated easily for small A and we obtain (Lf(0) = —2) 



Po(fc,0) = Po 



+ 



32tt 2 



(x + 4) 



k - 



TV 
6~J 



(6.5) 
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where we used the point k = T/6\ instead of k = A to start the evolution, as we have 
explained in the first paragraph of this section. We define the renormalized couplings of 
the theory in the limit k — > as [] 

Po =Po(0,0) 
Xr =A(0,0) 

x R =x(0,0). (6.6) 

and conclude that 

'°©^° + 3^ +4 K£) 2 - (6 - 7) 

At non-zero temperature, the evolution in the low temperature region (k > T/6\) is 
identical to the zero temperature case. In the threshold region (T/6\ > k > T/62), the 
form of ti t 2(w,T) is not given by a simple analytical expression. For small A we neglect 
the running of A, x in this region and find 



61 J 16tt 2 
1 

16^2 



» ' r Z 




72 



Ox 
T 

2 ' J \0i 



x ( — , t) =x ( — ) = x R , (6.8) 



where 

1= dyyhlo,-) (6.9) 
Jye 2 \ y) 

and we have made use of the fact that tf(w, T) depends on T only through the combination 
T/k. The integral I can be evaluated numerically. For 6\ = 0.15, 82 = 0.4 we find 
/ = 19.97. Eqs. ( |6.8|) set the initial values for the evolution in the high temperature 
region. 



7. The running in the high temperature region 

In the high temperature region (k < T/62) the functions L\ 2 {w 1 T) are given by the 
simplified expression 



Op . We can rewrite eqs. ( |6.1 



in terms of effective three- 



dimensional couplings 



PoiKT) 



_ Po(k,T) 
T 



7 In the case that Goldstone modes are present (as for x = 0) the couplings are defined at some 
appropriate non-zero k. The same applies for non-zero temperature. This does not affect our results for 
small A. For a detailed discussion see ref. 101. 
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\'{k,T) =\{k,T)T 
g'{k,T) =g{k,T)T 

^^W)- 1 - (7 - i) 



The flow equations then read 



d -^= -v 3 k{3Ll(2X'p' ) + (l + x)Ll(xX'p' )} (7.2) 
— = -v^fXf {9Ll(2X'p' ) + (1 + x) 2 ^(xAV )} (7.3) 

2 

% = v 3 k^^^-±{Ll(2X'p' ) - Ll(xX'p' )} 

+ v 3 A;- 1 xA / {9L3(2AV ) + (1 + x) 2 L 3 2 (xX' p' )} , (7.4) 

with t> 3 = 1/87T 2 . Comparison with eqs. ( |4.14| )- p.l6| ) shows that the above equations are 
exactly the ones of the three-dimensional theory at zero temperature. In order to make 
their fixed point structure more transparent we define the dimensionless parameters 

K(k,T) 

X(k,T) 



Po(k,T) 


Po(k,T) 


k 


kT 


X'(k,T) 


X{k,T)T 


k 


k 


g'(k,T) 


g{k,T)T 


k 


k 



9 (k,T)^-^ = ^-^. (7.5) 
In terms of these we obtain the scale invariant form of the evolution equations 

— k — v 3 {3Ll(2Xn) + (1 + x)L\ (xXk)} (7.6) 



dn 

~di ~ 
dX 

~dt ~ 

dx _ 6i + f 
dt V3 k 1 - | 



A - v 3 A 2 {9i^(2Afi;) + (1 + x) 2 L 3 2 (xXk)} (7.7) 
{lI(2Xk) - L'KxXk)} 



+ v 3 xX{9L 3 2 (2Xk) + (1 + x) 2 LI(xXk)). (7.8) 

No explicit dependence on the scale k appears on the r.h.s. 

The first of the above equations defines a critical surface K cr = K cr (X,x). It consists 
of the points K cr for which the solution of eqs. ( |7.6|) - (|7.8|) approaches, for large negative 
t, a scaling solution with k, X and x independent of t. (For a weakly first order transition 
the scaling holds only approximately.) Every point on the critical surface is unstable in 
the k direction. Trajectories which start at k > K cr continue towards the region of large 
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k, in such a way that po(k,T) = n(k,T)Tk reaches asymptotically a constant value for 
k — > 0. As a result the renormalized theory settles down in the phase with spontaneous 
symmetry breaking. If the evolution starts at k < n cr , the flows cross the surface k = 
at some finite k s . From this point on the system is in the symmetric regime. In order 
to continue the evolution, we define appropriate parameters according to eqs. (|3.1| ). The 
evolution equations read 

d 'TTl ^ — 

^- = v 3 k{A + x)X'Ll{m 2 ) (7.9) 
°±- = -vaJfc-^lO + 2x + x 2 ){\'] 2 L 3 2 (m 2 ) (7.10) 

(JjL 

dx — 

— = Vz k~\x + l)x{x - 2)X'L 3 2 {m 2 ), (7.11) 

with A' related to A by eqs. ( |7.1| ). Comparison with eqs. ( f4.1| ) - ( |4.3|) shows that the 
above equations are the ones for the three-dimensional theory in the symmetric regime. 
We define the dimensionless couplings 

rh 2 (k,T) 



m 2 (k,T) 
X(k,T) 



k 2 

\'{k,T) _ X{k,T)T 
k ~ k 



g(kiT) ^ iim ^jhiK. (7 . 12) 

In terms of these quantities the evolution equations in the symmetric regime read 

-2m 2 + v- i (A + x)\L\{m 2 ) (7.13) 



dt 



^ = -\-v 3 (10 + 2x + x 2 )X 2 Ll(m 2 ) (7.14) 

(JjL 

dx 

— = v 3 (x + l)x(x-2)\L 3 2 (m 2 ). (7.15) 

We start the evolution in this regime at k = k s with m 2 (k s ,T) = and X(k s ,T), x(k s ,T) 
taking their values at the end of the running in the spontaneously broken regime. For 
k — > the evolution is stopped by the mass terms in the threshold functions L\ 2 and 
the theory settles down in the symmetric phase. Obviously the critical temperature T cr 
is related to n cr (for given A(T cr /6 , 2), x(T cr /02)) by 

LrS n cr . (7.16) 



9 2 

On the critical surface there are two fixed points with at least one attractive direction 
for the flow towards the infrared (k — > 0): 
a) the Cubic fixed point located at 

n c = 5.674 x 10~ 2 A c = 8.747 x c = 2, (7.17) 
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b) and the Heisenberg fixed point located at 

k h = 4.486 x 1(T 2 \ H = 15.265 x H = 0. (7.18) 

Both fixed points are attractive in the A direction, but the Cubic fixed point is unstable 
in the x direction while the Heisenberg one is stable. For fixed x there is also the infrared 
unstable Gaussian fixed point. It is located at 

kg = -v 3 (x + 4)L?(0) = ^(x + 4) A G = 0. (7.19) 

These fixed points are the same as the ones observed in section 3. The only difference 
lies in the parametrization. In section 3 we used an expansion around the origin of the 
potential and the evolution equations in the symmetric regime. For this reason the fixed 
points appeared for negative values of the mass parameter (the curvature at the origin), 
indicating a minimum away from the origin. In this section we rely on a parametrization 
around the minimum of the potential, which results in increased quantitative accuracy 
for the truncation. 

The flows on the critical surface are qualitatively similar to those in fig. 1 for x > 0. 
There are two disconnected regions: 

a) 2 > x > 0: The trajectories flow away from the Cubic fixed point and towards the 
Heisenberg fixed point. This region corresponds to a second order phase transition. 

b) x > 2: The trajectories flow away from the Cubic fixed point and into a region of small A 
and large x. Similarly to our discussion in section 3, we expect X(k, T) to become negative 
at some finite k. This indicates that the minimum of the potential becomes unstable (it 
turns into a maximum). Our crude truncation is not sufficient for the investigation of 
this situation, since the higher derivatives of the potential are important. A detailed 
study is presented in ref. |33| for the four-dimensional theory at zero temperature. The 
dimensionality of the theory is not crucial for the qualitative behaviour in this region 
of parameter space. The first term in the r.h.s of eq. ( |7.7|) (which is present only for 
the effectively three-dimensional theory) is not important for large x. The second term, 
which drives the dynamics, is the same as for the four- dimensional theory (with the 
replacement of v^Lf by v^L\ generating only quantitative corrections.) In ref. |33| higher 
derivatives of the potential are taken into account. Also a parametrization is used which 
simultaneously follows the evolution of the potential at its minimum at non-zero po(k,T) 
and at the origin. This permits the study of the global properties of the potential. It is 
found that during the evolution in the region x > 2 a second minimum appears at the 
origin which subsequently becomes the absolute minimum of the potential. This results 
in a discontinuity in the order parameter and a first order phase transition. At some point 
in the evolution, A turns negative and the minimum at non-zero po disappears. From this 
point on the deeper minimum at zero is the only minimum. During the whole evolution 
the higher derivatives stay positive guaranteeing that the potential remains bounded. We 
cannot reproduce the above picture within the crude truncation of a quartic polynomial 
for the potential. Instead we shall give in section 11 an approximate solution of the 
evolution equation for the potential in this region, which will demonstrate the existence 
of the first order transition. 
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8. The critical temperature 



A quantity which can be easily calculated from the discussion in the last two sections is 
the critical temperature for the phase transitions. From eqs. ( |6.8|) and (|7.5|) we obtain 



^• r )^ + iMi- /|( ** +4) 



I -• V ' =6 2 \r 



2 



x[-,T)=x R . (8.1) 



2 

For small A the critical surface, which separates the symmetric phase from the phase 
with spontaneous symmetry breaking, goes through the Gaussian fixed point given by eq. 
( |7.19| ). The critical temperature can be computed as the temperature for which k (jj^, T 
coincides with the Gaussian fixed point kq- This gives 

Tl C 



po x R + 4 



c- 1 



in 2 



02 




independently of X R . Substitution of the values Q\ = 0.15 , 9 2 = 0.4, / = 19.97, which we 
computed in part I of this section, gives 

C = 23.89. (8.3) 

We should point out that the above value for the critical temperature is not strictly valid 
for x R > 2. In this region the first order phase transition occurs for k (j^,T\ slightly 
above the critical surface. As a result the transition takes place at a temperature slightly 
lower than the one given by eq. ( |8.2| ). In the language of the effective three-dimensional 
theory, the distance from the phase transition can be parametrized, for small X R , in terms 
of 5n cr = k T) — kg- We establish the connection between this quantity and the 
temperature as 

T \ / 1 1 



$K cr = K[—,Tj - k g = 6 2 p l — - — } . (8.4) 



" cr , 



The critical temperature can also be calculated in high temperature perturbation the- 
ory through the perturbative expansion of the effective potential || and its generalization 
to non-zero temperature The calculation is straightforward and we do not present 

the details here. When the leading term in the high temperature expansion of the one- 
loop contribution to the effective potential is retained, the critical temperature is found 
to be 

T 2 24 

" (8-5) 



- cr 



po x R + 4 
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This value is in excellent agreement with our result. The slight discrepancy is due to small 
deviations of the form of Lf(w, T) that we have used from the exact expression, and could 
probably be removed by using lower 8± and larger 62. It is well known that the perturbative 
expansion of the effective potential breaks down near the critical temperature for a second 
or weakly first order phase transition, due to infrared divergences @. The surprising 
accuracy of the perturbative estimate for the critical temperature is due to the fact that 
the infrared divergences appear at temperatures |T — T cr \/T cr = 0(\r). For sufficiently 
small \r the location of the transition can be accurately computed. However, naive 
perturbative predictions for the details of the transition (even its order) can be misleading. 
This is the case for our model, for which the next terms in the naive high temperature 
expansion of the perturbative result fail to reproduce even the correct qualitative picture 
f\ (apart from a small region in parameter space which will be discussed at the end of the 
next section). We emphasize that our approach may also be used for large values of Xr 
not so easily accessible to perturbation theory. In this case should be replaced by the 
relevant exact point on the critical surface k ct . 



9. Second order phase transition and the critical be- 
haviour 

Let us briefly summarize the main results of the previous sections. We have considered 
theories which at zero temperature are in the phase with spontaneous symmetry break- 
ing corresponding to the AX regime. They are defined in terms of the classical (bare) 
parameters at the ultraviolet cutoff A ^> T . The renormalized parameters are obtained 
by solving the evolution equations from k = A to k — 0. In the low temperature region 
(A > k > T/81) there is no difference between the zero and non-zero temperature case for 
the evolution of the parameters of the theory. The first temperature effects are observed 
in the threshold region (T/61 > k > T/82). For small \r the values of the running param- 
eters at the beginning of the evolution in the high temperature region (k = T/6 2 ) can be 
expressed in terms of the renormalized parameters of the zero temperature theory. The 
relation is given by eqs. (|6.8| ). In the high temperature region (k < T/62) the character 
of the evolution is effectively three-dimensional and is governed by the fixed points of 
the three-dimensional theory. These are most transparent in terms of the dimensionless 
parameters defined in eqs. ( |7.5| ). The evolution equations are given by eqs. ( |7.6| )-( [7T8D . 
These equations define a critical surface (A, x), which is unstable in the k direc- 

tion and separates the phase with spontaneous symmetry breaking from the symmetric 
one. The system ends up in either phase depending on the values of the running param- 
eters at k = T/82- For small A the critical surface goes through the Gaussian fixed point 
given by eq. ( |7.19| ). As a result, for small Xr the crucial quantity is the distance from the 

8 Perturbation theory for gap equations [^5|l may lead to more reliable results concerning the order of 
the transition but will fail for critical exponents unless the effectively three-dimensional running of the 
quartic scalar coupling is properly included 
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critical surface 5n cr = k (^,T^ — kq. For 5n cr > the theory ends up in the phase with 
spontaneous symmetry breaking. For 5n cr < it ends up in the symmetric one. There is 
a direct connection between 5n cr and the distance from the critical temperature, which is 
expressed in eq. (|8.4|) . 



In this section we discuss the region in parameter space < x < 2 where the phase 
transition is second order. Sufficiently near to the critical temperature the evolution 
is governed by the Heisenberg fixed point. In fig. 4 we plot the numerical solution 
of eqs. ( |7.6| )-( |7T8|) for the evolution in the high temperature region, for a theory with 
zero temperature renormalized parameters Xr = 0.01, xr = 1 and critical temperature 
T c 2 r /p = 4.78. We display two trajectories, which start slightly above and below the 
critical surface (and therefore correspond to temperatures slightly below and above the 
critical one). We observe that the system flows towards the Heisenberg fixed point which 
is attractive in both the A and x directions. It stays around this fixed point for several 
orders of magnitude in t and then deviates towards either the phase with spontaneous 
symmetry breaking or the symmetric one. During the "time" t = ln(fc/A) that the system 
stays close the fixed point it loses memory of the initial conditions of the evolution. Its 
dynamics is fixed solely by the fixed point, which has a purely three-dimensional character 
(as we demonstrated in section 7). As a result we expect that the behaviour of the theory 
near the critical temperature is independent of the details of the zero temperature theory. 
It must display universal critical behaviour characteristic of systems with Heisenberg fixed 
point. As long as n(k, T) stays almost constant around its fixed point value kr- and k — > 
as t — > — oo, the minimum of the effective average potential evolves towards zero according 
to 

p Q (k,T) = K H kT. (9.1) 

If the temperature is equal to the critical one, the system never leaves the fixed point and 
p (0, T cr ) = 0. If the temperature is slightly below T cr , K,(k, T) eventually runs away from 
the fixed point and diverges, so that p (k,T) = n(k,T)kT reaches a constant non-zero 
value as k — > 0. This value corresponds to the renormalized minimum of the effective 
potential at non-zero temperature and we denote it by 

Po(T)=p (0,T). (9.2) 

For a temperature slightly above T cr , n(k,T) (and therefore p (k,T)) runs to zero at a 
finite k s . From this point on the system is in the symmetric regime and the appropriate 
evolution equations are given by eqs. ( |7.13| )-( [f.l5| ). We start the evolution in this regime 



at k = k s with m 2 (k s ,T) = and X(k s ,T), x(k s ,T) taking their values at the end of the 
running in the spontaneously broken regime. For k — > the evolution is stopped by the 
mass terms in the threshold functions L\ 2 and the theory settles down in the symmetric 
phase. We define the renormalized mass in the symmetric phase as 

m 2 R (T) = m 2 (0,T). (9.3) 

We also define the renormalized couplings in both phases as 

X R (T) = X(0,T) 
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x R (T)=x(0,T). 



(9.4) 



It is important to point out that, while the system is staying close to the fixed point, the 
coupling X(k,T) evolves towards zero according to 

\(k,T) = \ H t (9.5) 

As a result \r(T) goes to zero as the critical temperature is approached. Its strong 
renormalization near T cr provides the resolution of the problem of infrared divergences. 
The ratio A#(T)T/t?1r(T) does not diverge near the critical temperature, in contrast to 
XrT '/mn(T). (Here Xr is the renormalized coupling of the zero-temperature theory, which 
is approximately equal to the bare one for small couplings.) We shall not elaborate on 
this point, but we refer the reader to ref. || for an extensive discussion. We also mention 
that as the temperature deviates from the critical one the system spends less "time" 
t = ln(A;/A) near the critical point. Its flow deviates from the those depicted in fig. 4 at 
earlier stages. As a result the universal behaviour ceases to dominate. 

The behaviour of the renormalized theory at various temperatures is shown in fig. 5 
for zero temperature parameters Xr = 0.01, x R — 1 and critical temperature T c 2 /p = 
4.78. We observe that po(T) moves continuously to zero, indicating a second order phase 
transition. The mass m R (T) is zero at T cr and positive for larger temperatures. The 
quartic coupling Xr(T) stays close to its zero temperature value for most temperatures, 
but is strongly renormalized towards zero near T cr . The ratio of couplings xr(T) again 
takes its zero temperature value, unless the temperature is sufficienlty close to T cr for 
the flow to reach the Heisenberg fixed point. The universal behaviour near T cr can be 
parametrized by critical exponents, which we define similarly to ref. H as 



Po(T) oc(T c 2 r - T 2 ) 2/3 
m\{T) cx(T 2 - T 2 r 
Ah(T) (x(T 2 - T 2 cr f 

x R (T) oc(T 2 - T 2 cr y. (9.6) 

More precisely, 2/3 is given by the derivative of lnp (^) with respect to ln(T c 2 — T 2 ), and 
similarly for the other parameters. The definiton of ( and fi applies only to the sym- 
metric phase. The exponent /3(T) is plotted in fig. 6 along with xr(T) for temperatures 
approaching T cr , for a theory with Xr = 0.2, xr = 1. It is apparent from the temperature 
dependence of xr that near T cr the Heisenberg fixed point becomes important. During 
its evolution the system stays long enough on the critical surface for this fixed point to 
generate universal critical behaviour. The exponent /3(T) approaches a temperature inde- 
pendent value which is independent of Xr and xr (as long as xr < 2) and characteristic 
of systems with Heisenberg critical behaviour. This value is 

p H = 0.32 (9.7) 
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in agreement with ref. |J. Two other exponents are fixed by the scaling laws and the 
finite value of the ratio \r{T)T / m R {T) . They are u H = ( H = 2j3 H . The above values 
for the exponents are in rough agreement with known values from three-dimensional field 
theory [^| |36fl . The agreement improves dramatically when less restrictive truncations 
are used for the study of the evolution equation for the potential, and wave function 
renormalization effects are taken into account flQfl . We have performed this more accurate 
calculation and obtained results which agree with the known values at the 4-5 % level. This 
work will be described in ref. |54) . Finally, the exponent /i also approaches asymptotically 
a constant value (c.f. eq. ( |7.8|) ) 



vh \ 6 
=0.64. 



L\{2\ h k h ) - L\(ti) 



+ A 



ii 



9L 3 2 (2\ h k h ) + Ll(0) 
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10. Tricritical point and crossover 

In the above discussion the Heisenberg fixed point was the only one which played any 
role. This was expected since the Cubic fixed point is repulsive in the x direction. Any 
flow which starts sufficiently far from it is further repelled and the system never feels its 
effect. However, it is possible that the values of the running parameters at the beginning 
of the evolution in the high temperature region are within the region of influence of the 
Cubic fixed point. An example is given in fig. 7, for a theory with A# = 0.01, xr slightly 
smaller than 2, and T 6 2 r /p = 3.98. As we have discussed in the introduction and section 
3, flows that start on the surface x = 2 in parameter space never move out of it. The 
flows depicted in fig. 7 start with x (jj>^) = 2 — 5x and fa < 1. For this reason, 
their deviation from the surface x = 2 is very slow. We display two trajectories which 
start a small distance 6n cr above and below the critical surface (and therefore correspond 
to temperatures slightly below and above the critical one, according to eq. ( |8.4| )). For 
\SkctI < 5i C 1 the flows stay on the critical surface and close to x = 2 for a large 
initial part of the evolution. During this "time" they approach the Cubic fixed point 
and stay near it. Finally x(k,T) starts growing and the system moves away from the 
repulsive (in the x direction) Cubic fixed point and towards the Heisenberg one. After it 
approaches this attractive (in the x direction) fixed point the evolution is similar to the 
one depicted in fig. 4. Systems which start with larger values of \^K cr \ behave similarly 
to fig. 7, but deviate from the critical surface at earlier stages of the evolution. As a 
result, they can feel the influence of both the Cubic and Heisenberg fixed point, or only 
the Cubic one, or they can deviate from the critical surface too soon for any universal 
behaviour to be induced. We calculate the renormalized parameters of the theory (at 
various temperatures) similarly to the previous subsection. Their behaviour as a function 
of temperature is analogous to that in fig. 5. The main difference concerns the small 
region around T cr . In this region the temperature dependence should reflect the influence 
of the two fixed points during the evolution. We first concentrate on values of 5n CT for 
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which the critical behaviour is dominated by the cubic fixed point. For this region we 
plot in fig. 8 the critical exponents corresponding to Po(T) and xr(T), which are defined 
according to 

p (T)(x(T c 2 r -T 2 )^ 
2 — xr{T) oc(T c 2 r — T 2 )~ Lp . (10.1) 

We observe that they reach constant values as the critical temperature is approached. 
The value for (3 should be characteristic of the Cubic fixed point. We find 

(3 C = 0.25 (10.2) 

and vq = (c — 2/3cs m agreement with the scaling laws and the finite value of the ratio 
\r(T)T /mji(T). We expect the Cubic fixed point to generate the universal behaviour 
characteristic of an Ising system. This is due to the fact that the theory decomposes into 
two disconnected ^-symmetric theories for x = 2 (see introduction and sections 3 and 
4). Indeed, the critical exponents that we have calculated are in exact agreement with 
the results of ref. for N = 1, which were obtained at the same level of the truncation 
scheme. Improved truncations result in values for the exponents which are in agreement 



with three-dimensional field theory |36| at the few percent level [33[. The exponent 
if = 0.16 is a typical example of a crossover exponent [^, p5|] . It is related to the growth of 
the unstable coupling at the Cubic fixed point, and therefore to the negative eigenvalue of 
the matrix which governs the evolution of small perturbations around the fixed point value 
of the parameters. We postpone a more detailed discussion of the crossover behaviour for 



a future publication p4 



The behaviour shown in fig. 8 changes if the critical temperature is further approached 
(extension of the graph to the right). Eventually the system moves away from the Cu- 
bic fixed point and the exponents (3, u, £ take values different from those typical of an 
Ising system. Also the temperature dependence of 2 — xr(T) cannot be described by a 



crossover exponent anymore and xr(T) will rather follow eqs. (|9.6|) , ( |9.8|) . We display 
this behaviour in fig. 9. The values of xr(T) give an indication of which fixed point in- 
fluences the system. It is clear that the Heisenberg fixed point takes over from the Cubic 
one very close to T cr . The temperature dependence of the exponent (3 is a characteristic 
example of a crossover curve. It demonstrates how the critical dynamics changes from 
Ising-like (for In ( t2 ~PA ~ -20) to Heisenberg-like (for In ( T2 ^ Qr ) ~ -50). A detailed 
discussion of this behaviour within more accurate truncation schemes will be given in ref. 

m. 



11. First order phase transition 

We turn now to the region x > 2 where we expect a first order phase transition, as 
we have explained in sections 3 and 7. Our truncation scheme is too crude to describe 
the behaviour of the potential in this region. We have approximated Uk(pi, P2, T) by a 
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second order polynomial in pi^. This permits the discussion of potentials with only one 
minimum. The study of first order transitions requires the use of improved truncations, 
where higher p-derivatives of Uk are taken into account and the possibility of two distinct 
minima is permitted. This has been done in ref. [33f for the zero temperature theory, 



and the existence of a first order transition has been established. We shall not repeat this 
calculation here. Instead we shall derive an explicit solution of the evolution equation in 
the region of large x, which will demonstrate the existence of first order transitions for 
the high temperature theory. 

In fig. 10 we plot the numerical solution of eqs. (|7.6[) - (|7.8|) in the high temperature 
region, for zero temperature renormalized parameters \r = 0.01, xr = 2.01, 3, 5. The 
temperature is very close to the critical one. We notice that for all three sets of parameters 
the evolution leads to a region of large x. In fig. 10 the curves for k and A are terminated 
when x = 30. We observe that the running parameters tend towards the same area of 
parameter space. More specifically, for x = 30 we find (very roughly) A ~ 3, k ~ 0.08. 
This convergence of flows was already apparent in fig. 1. The difference in the evolution 
lies in the "time" t = \n(k/A) that it takes for the various flows to reach the same region. 
The flows (a) and (b) are fast, while the trajectory (c) starts very close to the surface 
x = 2, is first attracted towards the Cubic fixed point, and finally deviates towards the 
region of large x. The Cubic fixed point separates the region x < 2, where we have 
observed second order phase transitions, from the region x > 2, for which we expect first 
order transitions. For this reason it is characterized as a tricritical point. (The Ising fixed 
point exhibits similar behaviour.) 

In the regions of large x we have g ^> A. As a result, the contribution of the 0i 
fluctuations to the evolution of t4(pi, 0, T) is suppressed as compared to the contribution 
of 02- Moreover, the increase of x in this region is mainly due to the fast decrease of 
A = Xk/T. In contrast, the coupling g evolves only slowly. The pi-dependent mass term 
for the 02 field is approximately given by gp\ (for P2 = and apart from a very small 
region around the origin) and has again a mild fc-dependence. Let us assume that for a 
given scale ko the solution of the truncated evolution equations depicted in fig. 10 gives 
a good approximation to the exact solution for the potential. (This means, in particular, 
that a two-minimum structure has not appeared yet at this scale for the true potential.) 
We denote the parameters of the theory at the scale k by k = n(k ,T), A = X(k ,T), 
9o — g(ko,T), xo = x(k ,T), and the mass term for the 02 field by goPi- Based on the 
remarks at the beginning of this paragraph we can obtain in the high temperature region 
an approximate solution of the evolution equation ( |5.3| ) for the potential on the pi-axis 
(p 2 = 0). By neglecting the first term in the r.h.s. of eq. (|5.3|) and the /c-dependence 



of U2 = Ip*-, the differential equation (5^3) is easily integrated. We obtain in the limit 
k — > (up to an irrelevant constant) 



U( Pl ,0,T) =f4 =0 (pi,0,T) 



1 - T r°° 
~A (pi - K k T) 2 - — / dxy/xhi 

2 87T Z JO 



Pkpjx) + g pi 
x + gopi 



(11.1) 
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The effective inverse propagator P(x) is given by eq. ( 2.11| ) and we have indicated that it 
must be evaluated for k = k . Together with the numerical solution of the flow equations 
near the critical surface for k > ko, which provides the "integration constants" Ao, (fa and 
KokoT, we expect the effective potential of eq. (|11.1| ) to be a very good approximation. 
(For a sufficiently small ratio of couplings Xr/qr we may identify ko with T/6 2 . This 
essentially reproduces the results of high temperature perturbation theory.) 

The effective potential of eq. ( |1 1 . 1|) describes indeed a first order phase transition. 
This can be most easily visualized if we approximate for the purpose of demonstration 



fco 



-X 



-k 2 



for x > kl 
for x < /cq- 
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One finds for the pi-derivative 
dU( Pl ,0,T) 



Ux(pi,0,T) 



dpi 



- K \ k T + X pi + -^-T I " dx^fx 



JO 



l :1 



1 



1 



x + g pi kl + g p! 



(11.3) 



Using a rescaled field variable 



this yields (with A = A T//c 
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For k < ka = the potential U(p\, 0,T) develops a minimum at the origin (pi = 0). 
For k only slightly below ka the origin is only a local minimum whereas the global 
minimum occurs at p\ ^ and the model is in the phase with spontaneous symmetry 
breaking. For sufficiently small k /ka, however, the absolute minimum is at the origin 
and the model is in the symmetric phase. (Note that | — v /parctan + §•[+■= is a 
positive function for all p.) There is a critical ratio k /ka (depending on the value of 
Ao(l + xo) 2 ) for which the minima at p\ — and p\ ^ are degenerate in depth, but 
they are still well separated from each other. Changing kq (which is a function of T) 
through this critical value leads to a first order phase transition with a jump in the order 
parameter. 

The necessity of a first order phase transition can also be seen by considering the 
pi-dependent quartic coupling Uu(pi) = 9 U ^j°' T ^ which obeys 



U n (pi 



\ 90 rp 

8n 2 



dx\fx 



{x + goPif 



(P ko (x) +g pi) 2 



;il.6) 
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By keeping only the most singular behaviour of the integral for p\ — > we obtain 

-2 T -3/2 1 

u 1 M = h + $ir-W- T ^=- (11 - 7) 

37T Z K i07T Jpx 

We have recovered the leading perturbative result for the behaviour of the quartic coupling 
near a first order phase transition. If the minimum of the potential pio(k) is sufficiently 
close to zero at the scale ko, the remaining evolution of pw(k) from ko to k = causes C/n 
to vanish at some scale k between and k . As a consequence, the minimum at p\o 7^ 
becomes a saddlepoint and disappears subsequently. Already before, a new minimum 
has been generated at the origin, which remains the only minimum in the subsequent 
evolution to k — 0. Since Uu(pi) is always negative for sufficiently small p\, the phase 
transition can never be second order and all values x > 2 must lead to a first order phase 
transition. 

Let us finally discuss a suitable choice of the scale ko from which on we can replace 
the numerical solution of the flow equations ( |7.6|) - ( |7.8|) by the approximate solution 
given by eq. ( |11.1| ). On one hand x must be sufficiently large in order to justify the 
neglection of the contribution of the ^-fluctuations in the approximate solution. On 
the other hand ko should be sufficiently high so that a second minimum at the origin 
has not yet been generated and the truncation of a polynomial around po is still valid. 
This requires that trajectories near the critical trajectory not end at k = too deeply 
in the symmetric phase. A realistic choice of ko should rather correspond at k = to 
the situation where two minima exist simultaneously. For the "quasicritical" trajectories 
depicted in fig. 10 a reasonable compromise for ko seems to be given by the value for 
which x reaches 30. (This corresponds to ka — 0.6.) The trajectories (a), (b) and (c) 
shown in fig. 10 correspond then to potentials U(p\,0,T) with two different minima, as 
can be seen from fig. 11. They are close to, but not equal, to the critical trajectories 
for which k would deviate from fig. 10 towards the end of the running, thus leading to a 
potential with two degenerate minima. 

In summary, we have established the occurence of a first order phase transition for 
x > 2. Moreover, we have reproduced the perturbative prediction for the form of the 
potential near the origin. We should emphasize, however, that the perturbative expression 
applies only to the integration of fluctuations from the scale k (at which x 3> 1) to zero. 
The flow from the region of x near 2 to the region where the perturbative expression 
becomes valid can be computed only through the use of evolution equations. The different 
flows correspond to first order transitions of varying strength. The discontinuity in the 
expectation value is of the same order as ko. Also the mass gap at the critical temperature 
is proportional to this scale. In consequence the discontinuities in p for the flows (a), (b), 
(c) in fig. 10 have a ratio of Ap a /Ap b /Ap c = 1/0.016/5.1 x 10 -9 . The last flow, which is 
remains in the vicinity of the tricritical point before deviating towards the region of large 
x, corresponds to an extremely weakly first order transition. 

Our results can easily be extended to the region x < 0. We have seen in the introduc- 
tion and section 4 that the AX regime (x > 0) and the M regime (x < 0) can be mapped 
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onto each other through a simple redefinition of the fields (see eqs. fl4.12|) , Q4.13| )). For 
this reason, the physical behaviour in the two regimes is the same. For example, the 
Cubic and Ising fixed points generate the same universal behaviour, characteristic of a 
^-symmetric scalar theory. Similarly, a first order phase transition occurs in the region 
as < — 1. We shall not repeat our discussion for as < 0. All our results can be extended to 
this region by the redefinition of fields and couplings according to eqs. ( |4.12|) , (|4.13| ). 

12. Conclusions 

We have used the formalism of the effective average action for the study of the high 
temperature phase transition for a theory of two real scalar fields Xi,2, with the symmetry 
(Xi -Xi,X2 <-> -X2,Xi X2), and quartic potential 

V(xu X2) = \m\x\ + xl) + \MxI + Xl? + \xMxl (12-1) 

The phase diagram of the theory is divided into four disconnected regions: x > 2, 
2 > as > 0, > x > — 1, x < — 1. Three fixed points with at least one infrared sta- 
ble direction exist on the surfaces separating these regions: The Heisenberg fixed point 
(as = 0) is attractive in the A and x directions, and corresponds to a theory whose sym- 
metry is increased to 0(2). The Cubic (as = 2) and the Ising (as = —1) fixed points 
are attractive in the A direction and repulsive in the as direction and correspond to two 
disconnected Z 2 (xi,2 *-* — Xi, 2) -symmetric theories, which are equivalent. The model has 
a second or first order phase transition, with critical temperature well approximated by 
the perturbative expression if A is small. 

Theories with classical parameters in the regions 2>as>0,0>x>— 1 have a 
second order phase transition. Very close to the critical temperature the behaviour of 
the system is universal. It is characterized by critical exponents, which are determined 
by the Heisenberg fixed point. For theories with classical parameters near the surfaces 
as = 2, x = — 1 the influence of the Cubic or Ising fixed point can be observed near - 
but not too close to - the critical temperature. This leads to a crossover phenomenon, 
characterized by a crossover exponent and a crossover curve, for temperatures approaching 
T cr . The universal behaviour is initially determined by the Cubic or Ising fixed point for 
small enough (T — T cr )/T cr . As the critical temperature is further approached the more 
attractive Heisenberg fixed point dominates. No part of this rich structure associated with 
the second order phase transition can be observed within perturbation theory. We should 
mention that for small values of A the region in temperature where these phenomena 
appear is rather narrow. This changes for larger A, where the critical behaviour extends 
over a larger temperature domain without changing the universal results. Even though 
we concentrated in the present paper on small values of A for the purpose of comparing 
with analytical results, our method applies equally well to large A. 

A first order phase transition is observed in the regions as > 2, x < — 1. Therefore, 
the Cubic and Ising fixed points are tricritical points separating regions of second and 
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first order transitions. The perturbative expression for the effective potential is a good 
approximation only for x ^> 2 and x ~ —2. All theories near the critical temperature 
with classical couplings x > 2 or x < -1 correspond to renormalized theories with x 2 
or x ~ —2 at scales of the order of the mass gap of the model. However, we distinguish 
two classes of theories: 

I) For classical parameters x > 2 or x ~ —2 one finds a strongly first order phase 
transition. Here the effects of quantum or thermal fluctuations are well approximated by 
the perturbative expression for the effective potential. 

II) For classical parameters x ~ 2orx ~ -1 we predict a very weakly first order transition. 
The use of the renormalization group is indispensable for the correct incorporation of the 
quantum or thermal effects which strongly renormalize the theory towards the regions 
x > 2 or x ~ —2. 

Our results are relevant for multi-Higgs-scalar extensions of the standard model [pH[] 
and multi-scalar models of inflation [plj] . They cast doubts on the general validity of 



perturbative predictions for the high temperature behaviour of these models even in the 
case of small scalar couplings. High temperature perturbation theory was found to give a 
reliable estimate for the effective potential only in limited regions of the parameter space. 
Our non-perturbative method works for arbitrary values of the couplings in eq. ( |12.1 



With straightforward modifications - inclusion of x 6 couplings and anomalous dimension - 
it gives quantitatively precise predictions for all temperatures and all regions in the phase 
diagram. As demonstrated earlier || |Kj] these modifications are sufficient for computing 
critical exponents with a few percent accuracy. The model is easily extended to the 
case where \i an d Xi are iV-component vectors with internal SO(N) symmetries. The 
high temperature phase transition in other two-scalar models with a different structure 
of the potential - as for example supersymmetric two-doublet models - can be treated in 
complete analogy with the present work. 
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Appendix A: The integrals L 



The integrals L^{w) defined in eq. (|4.4[ ) , with the choice of eqs. 
infrared regulator, read 



L d n (w) = -2nl dyy 



4+i exp(-y) 



[1 -exp(-y)Y 



1 - exp(-y) 



+ w 



-(n+l) 



(E5I3n for the 



where u> = w/k 2 . In the text we use or synonymously for the integrals 

( |A.1|) , even though they depend only on the dimensionless ratio w. The integrals obey 
the relations 



-nL d n+1 (w). 



(A.2) 



For w > — 1, the integrals |L^(w;)| are finite, [] monotonically decreasing functions of w 
We define 

L*(0) = -2l d n (A.3) 
and give l d for the first three values of n 



I d 



H= 2(i-H)rg-i) 

# = 3 (l - 2 3 "l + 3 2 "i) r - 2 ) • 
The asymptotic expansion of L^w) for large arguments w — > oo is 



(A.4) 



Li(w. 



+ n(n + l) 



'(H^HMH*"™- <a - 5) 



The pole structure of L^w) is not relevant for this work. For a discussion see ref. Q. 
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Figures 

Fig. 1 Flows on the (X,x) plane for the three-dimensional theory. The evolution is deter- 
mined by eqs. flOD , C|378|) with m 2 = 0. 

Fig. 2 The integrals L\{w), L\{w). 

Fig. 3 Lf(w,T) I L\{w) as a function of T/k, for various values of w/k 2 . 

Fig. 4 The evolution of k, A, a; in the high temperature region, for temperatures slightly 
above and below the critical one, and Xr = 0.01, xr = 1, T c 2 /p = 4.78. The system 
approaches the Heisenberg fixed point before deviating towards the symmetric phase 
or the phase with spontaneous symmetry breaking. 

Fig. 5 The phase transition for a theory with Xr = 0.01, xr = 1, T c 2 /p = 4.78. 

Fig. 6 The critical exponent /3(T) and the parameter xr(T) as the phase transition is 
approached, for a theory with Xr = 0.2, xr = 1. 

Fig. 7 The evolution of k, A, x in the high temperature region, for temperatures slightly 
above and below the critical one, and Xr = 0.01, xr slightly smaller than 2, 
T 2 r /po = 3.98. The system approaches first the Cubic and then the Heisenberg 
fixed point, before deviating towards the symmetric phase or the phase with spon- 
taneous symmetry breaking. 

Fig. 8 The critical exponents (3(T) and <p(T) as the phase transition is approached, for a 
theory with Xr = 0.5 and xr slightly smaller than 2. 

Fig. 9 The critical exponent /3(T) and the parameter xr(T) as the phase transition is 
approached, for a theory with A^ = 0.2, xr = 1.99. 

Fig. 10 The evolution of k, A, x in the high temperature region for temperatures close to 
the critical ones. The zero temperature parameters are: 

a) Xr = 0.2, Xr = 2.01, 

b) Xr = 0.2, Xr = 3, 

c) Xr = 0.2, Xr = 5. 
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Fig. 11 The effective potential U = 4n 2 U /k^T as a function of = y/2g pi/k , resulting 
from the integration of eq. (|11.5| ) with k = 0.08, Aq = 3 and xq = 30. 
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